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ABSTRACT 


Cerenkov radiation is calculated for electron beams 
which exceed the velocity of light in a nondispersive 
dielectric medium. The electron beam is assumed to be 
bunched as emitted from a travelling wave accelerator and 
the emission region is assumed to be finite. 

The direction of emission of the Cerenkov radiation is 
related to the index of refraction which in turn is related 
to atmospheric index of refraction. The calculated index of 
refraction in the atmosphere predicts changes to_ the 
Cerenkov radiation pattern and the emission threshold 


condition. 
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I. INTRODUCTION 


Cerenkov radiation is the radiation produced by a charge 
or group of charges, moving at speeds greater than light in 
a particular medium. It has been investigated, starting 
with the experiments of Cerenkov in 1934 and described in 
terms of a charged particle moving faster than the velocity 
of light in a medium by Frank and Tamm in 1937. A summary 
of work to 1958 is contained in the book by Jelley [Ref. 1]. 

The concept of Cerenkov radiation has been used 
occasionally as a detector of high speed charged particles. 
A detector using this concept will selectively detect only 
charged particles exceeding the velocity of light in the 
medium through which they are traveling. Other than this 
application to particle detection, it has been considered as 
a curiosity, and the phenomenon has been usually thought to 
be a weak source of broadband radiation. 

Until recently the study of Cerenkov radiation has been 
restricted almost exclusively to radition in the optical 
region, and most theoretical and experimental studies have 
been devoted to the properties of Cerenkov radition in the 
Optical band. 

The Cerenkov angle is defined by cos8 = c/v = 1/nA, 
where c is the velocity of light in the medium, v is the 
velocity of the charged particle, n is the index of 
refraction, and 4 is the usual notation for the ratio of v 
to cy, where c,is the velocity of light in vacumm. The ever 
refraction in atmosphere is used to describe the capability 
of the atmosphere to bend an electro magnetic wave toward 
the region of higher refractive index. When conditions 
correspond to the standard atmosphere, a wave initially 
horizontal does not travel along a straight path but moves 
on a path bent slightly downward, but still not equal to the 


curvature of the earth except in unusual circumstances. 
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The atmospheric index of refraction is a function of 
temperature, pressure, and humidity [Ref. 2]. 

The purpose of this thesis is to understand how the 
Cenrenkov radiation produced by a charged particle beam 
propagating in air is affected by the the changes of index 
of refraction of the atmosphere. The atmospheric density 
range considered is appropriate for a range in altitude from 
sea level to 10 Kn. 


I ee THEORY 
A. CHARACTERISTICS OF CERENKOV RADIATION 


A charged particle or a charge bunch in uniform motion 
in a straight line in free space does not radiate. The 
situation may be different, however, in a dielectric medium 
for which k>k,, €>1, and index of refraction of the medium, 


n, can be written as: 


(KE fe]* (eqn 2.1) 


n 


where k,is a permittivity of free space, k is a permittivity 
of medium, and €£ is a dielectric constant of the medium. 

If the charge is moving fast, that is, at a velocity 
greater than that of light in the medium, a dipole field 
results which will be apparent even at large distances from 
the track. Such a field will be momentarily set up by the 
charge bunch at each line element along the track, and in 
turn each molecule in the line element will then radiate a 
brief electromagnetic pulse. The radiation will be spread 
over a band of frequencies corresponding to the various 
Fourler componts of this pulse. 

In the general case, the radiated wavelets from all 
parts of the track interfere destructively so that, at a 
distant point, the resultant field intensity is still zero. 
However, if the velocity of the particle is greater then the 
phase velocity of the radiated wave in the medium, it is 
possible for the wavelets from all portions of the track to 
be in phase with one another so that, at a distant point of 


observation, there is now a resultant radiation field. 
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Pigure 2.1 Huygens Construction of Wavelet 


for Cerenkov Radiation. 


It will be understood from the Huygens construction shown in 
figure 2.1 that this radiation is only observed at a 
particular angle Q@ with respect to the track of the 
particle, namely that angle at which the wavelets from 
arbitary points such as pl, p2 and p3 on the track AB are 
coherent and combine to form a plane wave front, BC. This 
coherence takes place when the particle traverses AB in the 
Same time that the light wave travels from A to C. If 8 is 
defined as v/c, where v is the velocity of the particle, then 
in a time t the particle will travel a distance AB =Qq,t 
and the light will travel a distance AC = ( c¢,/n ). From 


this geometry we obtain 


iva) 


cos@, = 1/nB (eqn 2.2) 


which is known as the Cerenkov relation, where @® is the 


direction of emission of the Cerenkov radiation. 
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Figure 2.2 p vs Electron Energy. 
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It is seen that; for a medium of given refractive index, 
n, there is a threshold velocity fp min = 1l/n, below which no 
radiation takes place. Figure 2.2 shows a_ relationship 
between A (i.e., v/c,) and electron energy in Mev. Note 
aac, an the case of air, l/n = 0.999732, to achieve 
Cerenkov radition the electron energy needs to be greater 
then 22 Mev. At this critical energy the direction of 
radiation of coincides with that of the particle( i.e., 
along the path of the particle). 

The above discussion applies to the case where the path 
of the charge bunch is infinite. For a finite path, the 
phasing condition between the charge and the radiation is 
relaxed so that the sharp Cerenkov radition peak is 
broadened over a range of angles. The resulting radiation is 


given by 
W(v, n) # Vo7QRe (eqn 2.3) 


for a beam of periodic bunches of fundamental frequency ; 


and by 
E(v, n)dv = Q R@ dv (eqn 2.4) 


for a single charge bunch. Here W is radiated power in 


watts/sr at the integral harmonic frequency 
v = 4% (eqn 2.5) 


and E is the energy radiated per unit solid angle per unit 


frequency range(Joules/sr/Hz). In either case, Q is the 





constant 
Wexehe 
Q = (eqn 2.6) 
Bne 


1 


where A is the magnetic susuptibility of the medium) scm 
the velocity of the radiation in the medium, and gq is the 
charge in a single bunch. The radiation function R is the 


same in both cases 


R = kL sin® I(u) F(k) (eqn Zea 


where 
I(u) = (sin u)/u (eqn 2.8) 


where k is wave vector, L is the finite propagating distance 
of the electrons, F(K) is a dimensionless form factor, i.e., 
the Fourier transform of the charge bunch is qgqF(K), and the 
radiation parameter u is given by u = kL(cos@ - cos@ )/2. In 
the forward direction, the radiation function (and thus the 
radiation)is dominated by the diffraction function TG 
which has a maximum at u = O and adjacent zeroes at (+)z and 
(-)7% - In the variable © , the maximum corresponds to Sij=mes 


and the zeroes correspond to 


COSG, = (ng)7! + pol Cen) (eqn 2.9) 


cCOoSOD = (n8)7! - 771 ee (eqn 2.10) 


respectively, where 4 = L/x and A is the wavelength in the 
medium of the rediation. 

For a infinite path, 9% = 98, and the radiation peak is 
very sharp. For a finite path, 9, ana ¢j differ more and 
more as the path becomes shorter. The onset of radiation is 
defined as the condition when 9, just enters the physical 


range (cos@, = 1). 
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This situation may arise through various combinations of 
index of refraction, beam energy and beam length. 

The index of refraction of air is close to 1 and its 
variation with temperature or pressure is very small so that 
the effects of changes ang! the constant Q are 
inconsequential. Otherwise the index of refraction n enters 
the formalism only through the parameter u. Here, however 
the effects may not be completely neglible since the 
radiation patterns and the onset of radiation are sensitive 


moO U. 
Be INDEX OF REFRACTION FOR THE ATMOSPHERE 


To obtain a reasonable understanding of the index of 
refraction for the atmosphere, we must consider some basic 
physical principles. These considerations are necessary if 
one is to have an appreciation for the relative effects of 
pressure, temperature and water vapor. From previous 
discussion, the index of refraction was determined by the 
dielectric constant, ¢ , (n =f€ ). 


The atmosphere is a gas, and Debye derived an expression 





that relates the dielectric constant, &€ , of a gas to its 
basic properites. The equation is 
° 22 ul: 
e-l __ 4m A hmm ec 
eaes eu M&t+s i Gir! 
amount of permanent dipole 
material polarizability contribution 


The physical meaning of each term is indicated in the 
equation. The symbols are as follows, A = Avagadros number, 
f = density, K = Boltzmans constant, M = molecular weight, 
4 = permanent electric dipole moment, t = the relaxation 
time required for external field-induced orientation of the 


molecules to return to a random distribution after the field 
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is removed,W= 21uf where f is the frequency of the external 
field. 

One concludes from Debye's analysis that for external 
fields with freqencies less than 100 GHz/s, WU<< 1, with the 
result that eqn 2.11 is written as 


e-l 47 AO fam ne (eqn "2. iz 
E+ 3M IKE 


This expression contains two of the influences of an 
externally impressed field decribed previously. The 
polarizability is due to the distortion of the molecules by 
the impressed field. It is noted that both non-polar and 
polar molecules experience distortion. The permanent dipole 
contribution is due to the orientation effect on the polar 
molecules only. 


Since n =/€ and since n = 1 we approximate 


€-1 | avel _ 2(n-1) (eqn 2.13) 
“ee nae 3 


Uist na the above expression we can write the 
contributions to n-1 from non-polar (dry) and polar (Hae 


gases. We also use the equation of state / = P/RT. 


(eqn 2.14) 
eek 
(nl = >, AL _d 
lary "LR tee 
A u 
(n—1) oy ie eee (eqn 2.15) 
age MR T {w 3 kT? 


ie 


It is noted that e is the pressure associated with the 
water vapor and p, is the pressure associated with dry air. 
The total pressure, p , is their sum (p =e e ). When we 
add the contributions due to water vapor (polar) and dry air 


(non-polar) we get the index of refraction of moist air: 


Zee 
(n-1) =¢ (eqn ) 


1 P-e le i 
1 RS) AS ae US aS ae, 
Where C, =H /3K and C, = 2mA/R. Note that all parameters 
in Eqn (2.8) with the exception of p,T, and e can be treated 
as constants. ( Mg andedg are values averaged over the 
molecules of the atmosphere. ) 
If we combine the terms and assign values to the 


constants we obtain 


(n=1) = (77.6 = - 5.62 + 3.75 x 10° £) x 10° (eqn 2.17) 
T 


where p and e are in millibars and T is in degrees Kelvin. 
It is very convenient at this point to introduce a new 


parameter N which is the " Refractivity ", 


ie= ( n-1 ) x 10 (eqn 2.18) 


From the Ref.( 3 ), this equation becomes 


288p 6.7e (eqn 2.19) 


N = __ ee ee ee 
eee Ol 0.002661) GPUs ie) 


Where p and e are in millibars and t is in degrees Celsius 
mecunjeeieas the refractive index, defined as the ratio of 


velocity of propagation of the electromagnetic wave in a 


at] 


vacuummtomthatmam agus Since electromagnetic waves travel 
slightly more slowly in air than in a vacuum, the refractive 
index is always slightly greater than unity. At the earth's 
surface the numerical value of the refractive index is 
usually between 1.000250 and 1.000400. The value N is thus 
between 250 and 400. Refractivity can be expressed as a 
function of the atmosphere pressure, temperature, and 


humidity by the Equation (2.19). 
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Figure 2.3 Refractivity Profiles in the Standard Atmosphere 
from Ref. 2. 


In Figure (2.3), refractivity is shown as a function of 
altitude for surface refractivities of 240 and 400 N-units 
for standard atmosphere conditions. In the calculations to 
plot curves, layer thicknesses of constant refractivity of 


50 meters were taken from the ground level to an altitude of 
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30 Km. The region above that was considered to consist of 


unity refractive index [Ref. 2]. 


Me, 


III. CALCULATION OF REFRACTION INDEX IN THE ATMOSPHERE 
A. TEMPERATURE IN THE ATMOSPHERE 


The temperaturer of the surface air is well Known at 
many places and various heights, from sea level up to about 
6 kilometers. But the temperature records obtained by the 
aids of kites and balloons, both manned and free, show that 
the mountain air temperatures generally differ materially 
from the temperature of the free air at the same height and 
latitude. 

The average temperature of the surface decreases 
approximately at the rate of 1°C per each 180 meters, 200 
meters, and 250 meters increase of height on mountains, 
hills, and plateaus, respectively. In the free atmosphere, 
however, the result is quite different.Here the decrease of 
temperature with increase of altitude, except at very great 
heights is, roughly, the same at most parts of the world. 
Figure(3.1) from [Ref.13] shows at thetemperature decreases 


linearly between sea level and an altitude of 10 Km. 


B. PRESSURE IN THE ATMOSPHERE 


Atmospheric pressure (air pressure) is the force of 
hydrostatic pressure of air acting on a unit area. In 
immobile air it is equal to the buoyancy of air and is 
balanced at every point by the weight of a vertical column 
of air of unit cross section lying over the level under 


consideration. 
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Pigure 3. 1 Temperature Profiles in the Atmosphere 
from Ref.13. 


In practice, pressure is usually measured in terms of 
the height of a column of mercury in a barometer, expressed 
in millimeters ( mmHg ). But the height of the mercury 
column in a barometer, as is well Known, depends not only on 
pressure but also on the temperature of the mercury, as well 
as on the acceleration of gravity at the point of 
observation, which varies with elevation above sea level and 
latitude. The reading of a mercury barometer must therefore 
be reduced to a standard temperature (usually O° C), and a 
Standard altitude and latitude. Figure (3.2) from [Ref.13] 


72, \t 


shows the pressure decreases almost linearly between sea 
level and an altitude of 40 Km. 
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Figure 3.2 Pressure Profiles in the Atmosphere 
from Ref.13. 


C. WATER VAPOR IN THE ATMOSPHERE 

Water vapor enters the atmosphere because of the 
evaporation of water from the earth's surface and spreads 
through the atmosphere by mixing. At the same time the 
maximum content of water vapor is restricted by its 
saturation pressure; condensation of the vapor takes place 


when this point is reached. 


Zee 


The processes of evaporation and condensation are 
accompanied by the absorption or generation of heat, which 
is reflected in the thermal regime of the corresponding 
layers of the atmosphere. Water vapor also plays a large 
part in the processes of absorption and emission of 
radiation. All this explains the tremendous role of water 
vapor in all phenomena in the atmosphere.The mean vapor 
pressure e depends exponentially on the height above sea 


level. 


2 
10° 8 (eqn 3.1) 


Here e and e, are the mean vapor pressure at height z and 
at the surface respectively. The coefficient f is found from 
observational data. For the lower atmosphere, i 1s about 
5000m. From this it follows that e decreases by a factor of 
MOeat a Height of SKm, a factor of 2 at 1.5 - 2 Km, etc, 
i.e., The vertical decrease of the vapor pressure is much 
faster than that of the overall pressure of the atmosphere. 
The concentration of water vapor also decreases 
(exponentially) with altitude. Thus the vapor content falls 
off very rapidly with altitude and becomes negligible when a 
level of about 8 - 10 Km is reached; at great altitudes the 


mom Ls, as a rule, very dry( Figure 3.3 ). 


D. CALCULATION DATA 

Table I and II show the parameter for the US standard 
atmosphere [Ref.3] and the German atmosphere [Ref.4]. Using 
these data in equation 2.19, the microwave index of 
refraction can be calculated [Ref.5] as shown in Table III 
and IV. These results are then used to assess the changes in 
microwave Cerenkov radiation resulting from propagation of 


an electron beam in the atmosphere. 


23 


J, &e9 


S 





Ottis Weel 
PRESS eine 


ELlegure-3. 3 Water Vapor Profiles in the Atmosphere 
from Ref.4. 
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TABLE I 
U.S. STANDARD ATMOSPHERE 


Temp. Pressure (mmHg) 


C'Cin) vapor Saturatior 
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TABLE [I 
GERMAN ATMOSPHERE (MEAN 1900-1912) 


Temp. Pressure(mmHg) 


("Cims) vapor saturatior 
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TABLE III 


INDEX OF REFRACTION IN U.S. 


Alt. Temp. 


(°C ) 
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INDEX OF REFRACTION OF GERMAN ATMOSPHERE 


Alt. 
(Km) 


Temp. 
(°C ) 
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Ev OF CHAR NDEX © R RA ON_ON RENKO 


A. EMISSION THRESHOLD FOR COHERENT CERENKOV 
We define the onset of the emission of Cerenkov 


radiation to be the situation when 8, begins to enter the 


physical range [Ref.6]. Then setting @, = O in eqn. 2.10 
gives 
nB = n (n+1)7! cocie — Ie) 


A plot of eqn.4.1 is shown in Figure 4.1. 
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Figure 4.1 Threshold of Emission. 
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As beam length increases, the product ng first rises 
rapidly and then asymptotically approaches the value unity. 
Any value of nB above the curve gives Cerenkov radiation 
with a peak position dependent on the beam length. That is, 
for values of n 2 > 1, the Cerenkov angle @® is in the 
physical range. For values of nA& between the curve and 
unity, @ is non physical but radiation with a well defined 
peak is still produced. 

Using the usual relation between A and y. Equation 4.1 
can be put in terms off, the value of Y necessary for the 


onset of Cerenkov Radiation. 


1/2 (eqn 4.2) 


aaa 





Yen) =i 


This gives the required. ened) aoa. a E for the onset 
of emission in terms of the index of refraction n and the 
beam length. 

Limiting values of equation 4.2 can be obtained for very 
long and very short beam lengths. For infinite beam length, 


Ye(n=2) = [1 - (eqn 4.3) 


n 
which decreases and approaches the value of 1 as n 


increases. 


Table III and IV show that the entire range of index of 
refraction is between the 10 Km value n = 1.000092 and the 
surface water saturated value of n = 1.000500 for extreme 
condtions. In the following, the effects of variation of the 
index between these limits on the size and position of the 
first radiation lobe, the threshold electron energy for 


emission, and the radiation patterns, are assesed. 
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Figure 4.2 First Diffraction Lobe Limits 
for 20 MeV Electron Bunches. 


The bounds ® and © for the first diffraction lobe are 
given by eqn 2.9 and eqn 2.10 with the principal maximum 
lying between. Figure 4.2, 2.33 4.4, show the first 
diffraction lobe angular limits @ and 9, as a function of 
beam length parameter for 20 MeV, 100 MeV, and 1000 MeV, 
electron bunches. In these figures, the solid line 


represents calculations for the minium index of refraction 
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Figure 4.3 First Diffraction Lobe Limits 
for 100 MeV Electron Bunches. 


in the atmosphere and curve of large dashes marks the 
calculations for the maximum index of refraction in the 
atmosphere. 

In figure 4.3, 8, and 8 are given by the upper and lower 
curve of each pair. In both cases, 98, is large for short 
beam lengths, and @ is nonphysical. As the beam length 
increases, ®& falls and @ becomes physical and increases. 
Then in the limit of very long path, both angles converge 


to the value @ which is different for air of (avetermne 
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Figure 4.4 First Diffraction Lobe Limits 
for 1000 MeV Electron Bunches. 


index of refraction ( Q@= 1.81° the for upper curve and @ = 
0O.78° for the lower curve ). 

At an energy of 100 MeV, an electron is already moving 
very close to the speed of light as shown in Figure 2.2. 
Consquently, increasing the electron energy to higher values 
has only a very small effects. This is shown in Figure 4.4 
which differs only slightly from Figure 4.3. However, at 
energies lower than 100 MeV the electron velocity is a much 


stronger function of energy and at 20 MeV the electron 
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velocity is significantly smaller than ¢c, . The vetfecroma 
this change are shown is Figure 4.2 where the dashed curve 
( air saturated with vapor ) is similar to the dashed curves 
in Figure 4.3 and 4.4 except that the limiting value Gia 
is much smaller. However, The behavior of the curves for 
the 10 Km index of refraction ( n = 1.000092 ) are 
different. For these values of the index and electron 
energy, np < 1 so that @ is non physical and 8 does not 
exist for any value of the beam length. The behavior of Q@ 
at short beam length is the same as that for higher 
energies, but at longer beam lengths © falls precipitously 
to zero so that there is no Cerenkov radiation in the limits 
of very long paths. 

The value of }/ necessary for ®, to become physical is 
given by eqn 4.2, and is calculated for electrons in Figure 
4.5. Again, the solid curve is for the upper atmosphere 
value of the index and the dashed curve is for the saturated 
value. At short beam lengths the curves are identical, but 
at longer beam lengths the electrons require a higher 
velocity (through an increase in energy) in order for ® to 


become physical. 
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Figure 4.5 Threshold Electron Bunch Energies. 


The interpretation of the curves in Figure 4.5 is as 
follows; If a given physical situation is represented as a 
point with coordinates E and Mon the graph, when that point 
lies below the curve, Cerenkov radiation is not produced. 
Above the curve, but below the threshold energy for infinite 
path, is the transition region in which the radiation peak 
continuously increases in size but changes position only 
Sergh tly. Above the threshold energy for infinite path, 
Cerenkov radiation is always produced, with an associated 


narrowing of the first lobe as the path increases. 


33 


The indices of refraction are 1.000092 for upper solid 
line and 1.000500 for the lower large dashed line. Until a 
beam length of about 7 = 100, is reached, the minimum index 
of refraction and maximum index of refraction do not affect 
the threshold energy in the atmosphere. After that value,the 
minimum index of refraction slightly more affects to the 


energy. 


B. RADIATION PATTERNS 

Figure 4.6 through 4.12 show the effects of changes in 
the index of refraction on the radiation function R. These 
are all calculated for parameters of the NPS Linac: a 
fundamental frequency of 2.85 GHz and assuming a Gaussion 
charge distribution within a rigid electron bunch but wae 
variable energy. 

The beam length parameter W = LA enters into R 
directly and through u. Using the frequency of 2.85 GHz 
leads to a wavelength in the medium (air) of 10.49 cm for 
the fundamental. All of the beam lengensi=1 =) are calculated 
using this value of wavelength. 

Figure 4.6 shows the calculated power W(watts/steradian) 
radiated from a periodic beam of electron bunches traversing 
a path (beam length) of 100.0 cm for media of three 
different values of index; n = 1.000268 n = 1.000092 and n= 
1.000500. In this case, all three radiation patterns are 
identical and show a broadening of the radiation pattern as 
a result of diffraction due to the small size of the source. 

Figure 4.7 is calculated for 100 MeV electrons and a 
path length of L = 10,490 cm (|= 10 ). At thats longer path 
length the radiation peak is narrowed and heightened, and 
the maximum value occurs at a smaller angle. Differences 
are readily apparent between the pattern for the different 
indices. Here, the dotted line represents a nominal sea 
level dry air value ( n = 1.000268 ) and the dashed and 
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solid line represents the 10 Km and saturated value 
respectively. 

The upper bound of the main lobe is given by eqn 2.9 
Inspection of this relation shows that 96, as well as Q 
increases with increaSing index of refraction. This 1s shown 
clearly in Figure 4.7 where the first zero of the radiated 
power is least for the lowest index and then increases with 
increasing index. The calculated values for 9, are 2.66°, 
2.87°, and 3.12° corresponding to the respective values of 
1.000092, 1.000268, and 1.000500 for the index of 
refraction. As ® increases the position of the radiation 
peak also increases and is correspondingly heightened by the 
poem O factor in the radiation function R 

Large changes in the radiated power occur when the 
parameters change so as to pass the emission threshold 
boundary delineated in Figure 4.2, 4.3 and 4.4. At a beam 
length of L =10.494 x 10 cm Gin 10 ) the calculated 
threshold value for electron boundes is 21.676 MeV for 
nominal sea level dry air, n = 1.000268. Small departures 
from these values of the parameters can allow the threshold 
to be crossed with dramatic results. 

Figure 4.8 is shows the calculated radiated power as a 


function of angle for electron bunches having a beam length 


Seo) = 10 and an energy one percent below the above 
threshold of 21.676 MeV. The dotted line calculatedusing 
the nominal sea level value of the index ( n = 1.000268 ) 


shows sub-Cerenkov radiation in an oscillating spatial 
pattern but without the large principal radiation lobe. 
Increasing the index to n = 1.000282 ( appropriate for sea 
level with 2.49 mb water vapor at US standard atmosphere) 
clearly shows the onset of Cerenkov radiation and the 
increaseing size of the main lobe. Further increase of the 
index drives the curve higher and to larger angles as a 


increases. 
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Figure 4.9 shows the radiated power as a function of 
angle for the same beam length and energy but for the 
extrems saturated water value of the index ( n = 1.000500 ). 
Here the peak power is 15 x 10 W/sr and it impossible to 
show the curves from the other indices on the same scale. 

Figure 4.10 is plotted to an intermediate scale and does 
show both the n = 1.000282 and the n = 1.000500 curves 
although the peak for the latter is far off scale. In Fig. 
4.9 and Fig. 4.10, the type of line associated with a 
particular index of refraction is different (f£0r alvgea 
altitude reasons) from the previous figures. When the 
radiation parameters are far from the threshold values, 
changes in the index are less drastic. 

Figure 4.11 shows the effect of varying the index of 
refraction on 100 MeV electron bunches with a beam length of 
m = 100. For this rather short beam the differences in 
radiated power are small but readily apperent. 

The position and size of the peaks corresponding to the 
different indices are in the same order as Figure 4.7, 4.8, 


4.9 and 4.10 for the same reasons. 
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Figure 4.6 Power Radiated vs Angle(l). 


Power radiated per unit solid angle (W/sr) as a value of 
angle for 100 MeV electron bunches and the beam lengths for 
100 cm. The function represents three index of refraction n 
= 1.000092, n = 1.000268, n = 1.000500, i.e; all of which 
are identical. 
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Figure 4.7 Power Radiated vs Angle(2). 


Power radiated per unit solid angle as a value of angle for 


100 Mev electron bunches and a 10490 cm beam length ( 
), n = 1.000500 for dashed lines; on 
1.000092 for solid line. 


= 10 


= 1.000268 for dota 
n = 
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Figure 4.8 Power Radiated vs Angle(3). 


Power radiated per unit solid angle as a value of angle for 
21.459 MeV electron bunches and a beam length of 10.494 x 


meorem ( = 10>), n = 1.000268 for dot line, 


n = 1.000282 
bor Gashea dot line. 


The energy 21.495 MeV is one persent 
below the threshold for n = 1.000268. 
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Figure 4.9 Power Radiated vs Angle(4). 


Power radiated per unit solid angle as a value of angle for 
21.459 MeV, n = 1.000500. The energy and beam length values 


are the same as for Fig. 4.8. 
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Figure 4.10 Power Radiated vs Angle(5). 


Power radiated per unit solid angle as a value of the angle 


for the same condition of Figure 4.9. This Figure shows the 


radiation pattern for n = 1.000282 (dashed dot line) and n = 
1.000500 (solid line). 
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Figure 4.11 Power Radiated vs Angle(6). 


Power radiated per unit solid angle as a function of 
altitude for L = 1049.0 cm and E = 100 Mev. The nominal dry 
sea level for dotted line (n = 1.000273), 10 Km height for 
solid line (n = 1.000092), extreme water saturated sea level 
for dashed line (n = 1.000500). 


42 


V. RESULT AND DISCUSSION 


The main points of this paper are to asess, in the GHz 
frequency region of the electromagnetic spectrum, the 
effects of height and atmospheric water vapor on the index 
of refraction, then to determine how this index of 
refraction affects the radiation pattern and the emission 
threshold of Cerenkov radiation. 

Caculation of the index of refraction for the atmosphere 
up to 10 Km as shown in Table III and Table IV are new 
results. The table values vary between 1.000092 and 
1.000500. These two values were used in calculating the 
thresholds, radiation pattern, and angular limits of the 
radiation pattern. 

The principal results of this study are; 

1 For high energy electron beams the principal 
radiation lobe is concentrated in a narrow range about the 
angle 8. as the beam length increases. The particular value 
of Q@ depends on the value of the index of refraction with 
larger indices giving larger values of &. 

2. For low energy beams, the lobe may or may not narrow 
depending on the value of the index. 

3. The threshold energy value for the onset of Cerenkov 
radiation at long beam lengths is strongly affected by 
changes in the index of refraction. 

For a specific example, using Fig. 4.5, the variation of 
index of an/n = 0.000408 gives an approximate doubling of 
the threshold energy. 

4. The radiation patterns for long beam lengths are 
sensitive to changes in the index of refraction. This 1s 
true for beam energies for above threshold (Fig. 4.7) and 
for beam energies near threshold (Fig. 4.8). 

These results lead to the following conclusion about an 


electron beam propagating in the atmosphere. 
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As an electron beam propagates upward in the atmosphere, 
it loss energy and sees a diminishing index of refraction. 
If the beam radiates over its entire length the effect will 
be to narrow the radiation cone and diminish the radiation 
as the diminishing beam propagates upward. If the beam is 
prevented from radiating by a plasma shield along most of 
its length, then only the short upper part of the beam would 
be expected to radiate. In this case, diffraction broadening 
could be severe, expecially for the lower frequency 


harmonics of the radiation. 
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